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1 Introduction 



In the standard calculations of ab oscillations in the medium [1,2] the one-particle (potential) 
model is used. The interaction of particles a and b with the matter is described by the potentials 
U a £. Im£4 is responsible for loss of 6-particle intensity. In particular, this model is used for the 
ab transitions in a medium followed by absorption of 6-particle: 

(a — medium) —)■(& — medium) — > f. (1) 

Here {b — medium) — > f represents the 6-particle absorption or decay. It is common knowledge 
that the models with non-Hermitian operators have very restricted area of applicability. In [3-5] 
this problem has been studied by the example of nn transitions in a medium [6-14] followed by 
annihilation 

(n — medium) — > (n — medium) — > M, (2) 

where M are the annihilation mesons. The reason is that for the process (2) the absorption 
in the final state (annihilation) is extremely strong. This is simplest 2-step process involving 
intermediate-state interaction and absorption in the final state since the nn transition vertex 
corresponds to 2-tail diagram. It was shown [3-5] that one-particle model mentioned above does 
not describe the process (2) as well as the total nn transition probability, because unitarity 
condition (optical theorem) is used for the essentially non-unitary ^-matrix. The interaction 
Hamiltonian contains the antineutron optical potential U n and lmU n plays a crucial role. The 
S'-matrix should be unitary. Since the annihilation is the main effect which defines the speed of 
process (2) in the medium, the potential model should be rejected. (Potential model describes 
only ab transitions with 6-particle in the final state [3] when unitarity condition is not used.) 

In [4,5] the field-theoretical approach with final time interval has been proposed. However, 
the model with bare propagator (see Fig. la), that is, the case with antineutron self-energy 
£ = has been considered only. In this paper the model with dressed propagator (see Fig. lb) 
is studied as well. It is shown that the result is very sensitive to the E. The uncertainty in the 
value of the lower limit on the free-space nn oscillation time r is conditioned by the uncertainty 
in the value of E. 

The principal object of this paper is to study the role of intermediate state interaction, 
namely, the role of particle self-energy induced by intermediate state interaction. The conclusion 
obtained for process (2) is generalized to process (1). For example, the possible AA transitions 
[15,16] in a medium followed by A decay: 

(A — medium) — > (A — medium) — > (n°n — medium). (3) 

2 Calculations 

The qualitative picture of the process (2) is as follows. The free-space nn transition comes from 
the exchange of Higgs bosons with the mass m H > 10 5 GeV [7] and so the subprocess of nn 
conversion is scarcely affected by a medium effects. From the dynamical point of view this is 
a momentary process: r c ~ l/rrin < 10~ 29 s. The antineutron annihilates in a time r a ~ 1/T, 
where T is the annihilation width of n in the medium. We deal with two-step process with the 
characteristic time r c h ~ r a . 

We consider the process (2). The neutron wave function is 

n(x) = Vt~ l l 2 exp(— ipx). (4) 
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Here p = (e, p) is the neutron 4-momentum; e = p 2 /2m + U n , where U n is the neutron potential. 
The interaction Hamiltonian has the form 

Hj = H nn + H, (5) 

H nii (t) = J d 3 x(e nn ty n (x)y n (x) + H.c), (6) 

H(t) = J d 3 xT-L(x). Here H is the Hamiltonian of the n-medium interaction, H nn is the Hamil- 
tonian of nn conversion [10], e nfl is a small parameter with e nfl = l/r nfi , where r nfi is the 
free-space nn oscillation time (r nn > r by definition); m n = m n = m. 

In the lowest order in H nn the amplitude of process (2) is uniquely determined by the 
Hamiltonian (5): 

M = e nn G M a , (7) 

°° = e fl -p2/2m-tf B + »0' (8) 

Pn = p, e n = e. Here Go is the antineutron propagator. The corresponding diagram is shown 
in Fig. la. The amplitude of antineutron annihilation in the medium M a is given by 

</0| TexpH / dxH(x)) - 1 |0n p >= A^(2vr) 4 5 4 (p / - Pi )M a . (9) 

Here | 0n p > is the state of the medium containing the n with the 4-momentum p — (e, p); <f\ 
denotes the annihilation products, N includes the normalization factors of the wave functions. 
The antineutron annihilation width T is expressed through M a : 

d<$>\M a \ 2 . (10) 

Since M a contains all the n-medium interactions followed by annihilation including an- 
tineutron rescattering in the initial state, the antineutron propagator Go is bare. Once the 
antineutron annihilation amplitude is defined by (9), the expression for the process amplitude 
(7) rigorously follows from (5). For the time being we do not go into the singularity Go ~ 1/0. 
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Figure 1: a nn transition in the medium followed by annihilation, b Same as a but the 
antineutron propagator is dressed (see text) 
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One can construct the model with the dressed propagator. For the density of the Hamilto- 
nian % we consider the model 



H(t) = J d 3 xU(x) = H a (t) + V, (11) 

where T-L a is the effective annihilation Hamiltonian in the second quantization representation, 
V is the residual scalar field. We include the scalar field V in the antineutron Green function 

Gi = Go + Go V Go+ ... = wg i- T7 = -I = -i. (12) 

In line with physical meaning of Eq. (12) (power series in self-energy) we denoted V — E, 
where E is the antineutron self-energy. The process amplitude is (see Fig. lb) 

M = e nn G d M b , (13) 

GdM b = G M a . The block M b = M b (H a , V) should describe the antineutron annihilation. 
Both Gd and M b depend on V. Because of this the block (amplitude) M b cannot be naturally 
obtained from the formal expansion of the T-operator 

TexpH / dx(H a + V^n)), (14) 

as is the case for M a (see (9)). An additional restrictions on the expansion of (14) are required. 

The models shown in Figs, la and lb we denote as the models a and b, respectively. If 
E — > 0, the model b goes into model a. In this sense the model a is the limiting case of the 
model b. We consider the model b. For the process width one obtains 

2 2 

T b = N 1 Jd$\M \ 2 = |f N, J 'd$\M b \ 2 = |f V, (15) 

T' = N ± J d& \ M b | 2 , (16) 

where T' is the annihilation width of n calculated through the M b (and not M a ). The normal- 
ization multiplier is the same for T b and V. We recall the amplitude M b is unknown. (The 
antineutron annihilation width T is expressed through the M a .) For the estimation we put 

M b = M a , V = T. (17) 

This is an uncontrollable approximation. Now the process width is 

r 



b T 2 V2' 

T nn n 

The time-dependence is determined by the exponential decay law: 



(18) 



W b (t) = 1 - e- r >* « |f F't = -£—t. (19) 

Equations (18) and (19) illustrate the result sensitivity to the value of parameter E. The 
physical meaning of T b and W b is standard: there is one-particle metastable state described by 
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wave function (4). The width and probability of two-step decay (2) are given by (18) and (19), 
respectively. 

Let r b and T nn are the lower limit on the free-space nn oscillation time calculated by means 
of model b and the experimental bound on the neutron lifetime in a nucleus [13], respectively. 
From the condition Wb(T nn ) < 1 we have 

Tnn >T b = ^^TT^n- (20) 

. For estimation we take V = 100 MeV and E = 10 MeV. Using the value T n , n > 1.77 • 10 32 yr 
obtained by Super-Kamiokande collaboration [13], one obtains 

r b = 1.2 • 10 9 s. (21) 

This value exceeds the restriction given by the Grenoble reactor experiment [8] by a factor of 
14 and the lower limit given by potential model [9-13] by a factor of 5. 

At the point E = Eqs. (19) and (20) are inapplicable. In this connection we return to the 
model a. The amplitude (7) diverges 

M = e nn G M a ~^ (22) 

Go ~ 1/0. These are infrared singularities conditioned by zero momentum transfer in the nn 
transition vertex. (In the model b the effective momentum transfer q = V = E takes place.) 
For solving the problem the field-theoretical approach with finite time interval [14] is used. It is 
infrared free. In the case of potential model, the approach with finite time interval reproduces 
all the results on the particle oscillations obtained by means of potential model (see sect. 5.2 of 
ref. [4]). This is the test of above-mentioned approach. However, our purpose is to describe the 
process (2) by means of Hermitian Hamiltonian because the potential model does not describe 
the absorption. 

For the model a the process (1) probability was found to be [4,5] 

W a {t) » W f (t) = e 2 nn t 2 , Tt » 1, (23) 

where Wf is the free-space nn transition probability. Owing to annihilation channel, W a is 
practically equal to the free-space nn transition probability. If t — > oo, Eq. (23) diverges just 
as the modulus (22) squared does. If E — > 0, Eq. (18) diverges quadratically as well. Equation 
(23) does not contain density-dependence (r— dependence) since it corresponds to the limiting 
case Ft ^> 1 which is realized for the nn transition in nuclei [4,5]. 

The explanation of the independence is simple. The process shown in Fig. la represents 
two consecutive subprocesses. The speed and probability of the whole process are defined by 
those of the slower subprocess. If 1/r «C t, the annihilation can be considered instantaneous. 
So, the probability of process (2) is defined by the speed of the nn transition: W a ~ Wf ~ t 2 . 

Distribution (23) leads to very strong restriction on the free-space nn oscillation time: 

T a = 10 16 yr. (24) 

Here r a is the lower limit on the free-space nn oscillation time calculated by means of model a. 
For the process (1) instead of (19) we have 

W &) « zra*- ( 25 ) 

T ab^b 

Here r ao is the free-space ab oscillation time, is the 6-particle self-energy, T is the absorption 
(decay) width of 6-particle in the case of absorption (decay). 
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3 Discussian and conclusion 



In this paper we don't touch the application of above-considered approach to the specific ex- 
periments on oscillations of ultra cold neutrons in storage vessels [17] and neutron disapearance 
into another braneworld [18]. Also we don't consider the calculations based on the diagram 
technique for direct reactions [19,20] since they are inapplicable to the problem under study 
[21]. (For example, the result T ~ l/E? t [19,20] (E n is the neutron binding energy) diverges 
when E n — > 0.) Our goal is to study the role of intermediate state interaction. 

As is seen from Eqs. (19) and (25), the self-energy plays the same part as the energy gap 
Un — U n in the potential model: it tends to suppress the processes (1) and (2). If E — > 0, 
Wb rises quadrat ically. So r h and r a can be considered as the estimations from below and 
above, respectively. In both models the antineutron propagators don't contain the annihilation 
loops since the annihilation is taken into account in the amplitudes M a and Mb, the interaction 
Hamiltonians Hj and unperturbed Hamiltonians are the same. The sole physical distinction 
between models a and b is the non-zero antineutron self-energy in the model b which is con- 
ditioned by intermediate state interaction. However, it leads to the fundamentally different 
results (see (21) and (24)). 

There is no rigorous recipes for calculation of £. One can adduce a great deal of arguments 
in support of the value £ = [4, 22]. In particular, in contrast to the model b, the model a 
describes the following limiting cases [22]: 

a) In the low-density approximation the model a reproduces the branching ratio of channels 
of free-space nN interaction: 

£ = -■ (26) 

Here a a and a t are the cross section of free-space nN annihilation and total cross section of 
free-space nN interaction, respectively; T a and T t are the width of process (2) and total width 
of nn transition in a medium. (I\ includes the channels with antineutron and annihilation 
mesons in the final state.) 

b) The model a describes the process (2) on a free nucleon: 

n + N -> n + iV -> M. (27) 

However, the small change of antineutron self-energy E = — > E = V ^ 0, or, similarly, effctive 
momentum transfer in the nn transition vertex converts the model a to the model b: W a — > Wb 
with Wb <C W a . This is because the process amplitude is in the peculiar point (see (7) and (8)) 
owing to zero momentum transfer. If £ changes in the limits 10 MeV > E > MeV, the lower 
limit on the free-space nn oscillation time is in the range 10 16 yr > r > 1.2 • 10 9 s. 

Finally, the problem is extremely sensitive to the E. This conclusion is also true for the 
process (1) since for this processes the E-dependence of the result is the same (see (25)). Further 
investigations are desirable. The self-energy of intermediate particle should be studied first. 
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